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Abstract. We present a self-contained analysis of the stationary radiative transfer equation 
in weighted L'' spaces. The use of weighted spaces allows us to derive uniform a-priori 
estimates for 1 < p < oo under minimal assumptions on the parameters. By constructing an 
explicit example, we show that our estimates are sharp and cannot be improved in general. 
Better estimates are however derived under additional assumptions on the parameters. We 
also present estimates for derivatives and traces of the solution and formulate a natural 
energy space, for which the data-to-solution map becomes an isomorphism. As a side result, 
we are able to prove uniform convergence of the source iteration for all 1 < p < oo without 
the assumption of positive absorption that is frequently used in the literature. 
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1. Introduction 

The equilibrium distribution of mutually non interacting particles propagating through a scattering 
medium is described by the stationary radiative transfer equation [11, 12] 



V • V0(r, v) + cr(r, v)0(r, v) = / fc(r, v', v)(/)(r, v') dv' + /(r, v) 

Jv 

where n(r) • v < 0. 



'V 
5(r,v) 



(1) 
(2) 



This system arises as a basic model, e.g., in radiation hydrodynamics [28], in reactor physics [11], in 
astrophysics [26], in climatology [16], or in optical tomography [3]. Depending on the application, the 
function (j), is called angular flux or specific intensity. The unit vector v = v/[v| denotes the direction 
and V G V is the velocity of propagation. The V operator only involves derivatives with respect to the 
spatial variable r G TZ. Particles enter the system through interior sources / or via a fiux g over the 
infiow boundary r_, where 



± •= 



{(r, v) edTZxV : ±n(r) • v > 0} 



and n(r) denotes the unit outward normal for r S dTZ. For physical reasons, the total cross-section a 
and the scattering kernel k are non-negative functions of their arguments. Let us further define the 
scattering cross-sections 

(7s(r,v) = / fc(r,v,v')dv' and a'g{r,v) — / fc(r, v', v) dv', 
Jv Jv 

which describe scattering from or into direction v. The difference a ~ a^ is the absorption parameter. 
In the following we shortly review some of the basic solvability results for (l)-(2) and recall the 



Date: April 25, 2013. 

Department of Mathematics, Numerical Analysis and Scientific Computing, Technische Universitat Darm- 
stadt, Dolivostr. 15, D-64293 Darmstadt, Germany. 



2 HERBERT EGGER AND MATTHIAS SCHLOTTBOM 

conditions they are based on: Solvability in L^ and L°° has been established in [10] under the sub- 
criticality conditions 

cr»(r, v) ,, , cr'(r, v) 

cr(r,v) cr(r,v) 

These imply that the scattering operator is a small perturbation of the differential operator on the left 
hand side of (1) and contraction arguments apply. Corresponding results in L^ for 1 < p < oo can be 
found in [2, 14]. A rather complete LP theory has been developed under similar conditions in [15] in 
the framework of semi-group theory. Note that a-priori estimates for the solution derived under these 
conditions typically degenerate when :^ — > 0; see also Theorem 1.2 below. In [32, 25], solvability in L^ 
was established provided that 

(T — (Ts > 0, cr — cr(. > 0, and cr > 0. 

Also some additional assumptions on the set V of velocities are required. Existence results in L^ were 
developed under these conditions in [2, 23, 17] by variational arguments. Note that the assumption 
cr > excludes the presence of void regions and that the a-priori estimates again degenerate when a — > 
0. Based on monotonicity arguments, existence of solutions in L^ was established in [27, 19], without 
the strict positivity assumption on a. A similar result was obtained in [4] under some restrictions on 
the set of velocities V. The existence in L^ for 1 < p < oo was actually announced in [15] but without 
proof. For velocities with uniform speed |v|, solvability in L^ was established without lower bounds on 
a in [18]. While the previous results are based on some sort of contraction principle, it is possible to 
obtain existence of solutions also via compactness arguments and Riesz-Schauder or analytic Fredholm 
theory [32, 31]. These results however do not lead to computable a-priori bounds. Let us finally also 
refer to [10, 15, 29] and [32, 5, 33, 6, 24, 9] for analysis of time dependent problems, results of spectral 
theory, and further references. 

The aim of this manuscript is to unify and generalize previous solvability results, to relax the conditions 
on the parameters, and to sharpen the a-priori estimates. We will present a self-contained L^ theory 
for stationary radiative transfer under the following assumptions: 

(Al) Let V C M'' be open and 7?. C K'^ be a bounded Lipschitz domain. 

(A2) cr : 7?. X V —> M is non- negative and a£ G L°°{TZ x V). Here i{r,v) denotes the length of hue 

segment through r in direction v completely contained in TZ; cf. e.g. [13, 2]. 
(A3) k : TZ X V X V ^ M. IS non-negative and measurable and 

(7 — (Tg > and cr — a'g > 0. 

All our arguments apply almost verbatim to more general velocity spaces equipped with some positive 
(T-finite Radon measure fi with /x({0}) = 0; see e.g. [33, 6, 30]. We use the assumption of an open set 
equipped with the Lebesgue measure mainly for ease of notation. The first of our two main results is 

Theorem 1.1. Let (Al)-(A3) hold and let Cp ^ -\\<Tsi\\L-^ + ^^^||cr^^||L~ < oo. Then for all 
1 < p < oo and all admissible data f , g, the radiative transfer problem (l)-(2) admits a unique 
solution (f) that satisfies 

l|ri0||iP(7JxV) < e^''(lK'"^/l|LP(KxV) + ll5l|LP(r_;|n.v|))- (3) 

Here ||u||ip(£).^) = ||w^/Pu||lp(u) = (/^ w|u|Pd(r, v)) ^ is the norm of a weighted space L^ space. 
Almost all solvability results mentioned earlier can be obtained easily from this theorem as special 
cases. It is the use of weighted norms that allows us to derive a-priori bounds which are uniform for 
all 1 < p < cxD. We will show in Section 6 that the a-priori estimate (3) is sharp and state estimates 
for the directional derivatives v • V0 and the traces 0|r+ in Section 8. 

Like many of the previous results, the proof of Theorem 1.1 is based on a fixed-point argument. We 
will establish the contraction property with a factor 1 — tt, where tt = 1 — exp {—Cp) can be interpreted 
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as the probability that particles leave the domain via the outflow boundary F^. As expected, tt gets 
smaller when increasing the size of the domain or the scattering cross-section. Particles may escape 
the system also by absorption. This case is covered by 

Theorem 1.2. Let the conditions of Theorem 1.1 hold, let cr > 0, and assume that for some v > Q 
we have \\(Js/(^\\L"=-(nxV) < 1 - ^^ and ||crs/cr||L^(7jxV) < 1 - t'- Then 

h'^^LPin^V) < t^"^lk'^"VllLP(7?,xV) +i^"'^ll3lUp(r_;|v.n|)- (4) 

Now u plays the role of a probability that particles get absorbed when interacting with the medium. 
Similar results can be found in [25] and [32] . The use of weighted norms is again essential to obtain 
uniform estimates for all 1 < p < oo, and it allows us to obtain also simple bounds for the directional 
derivatives and traces which will be stated in Section 8. In contrast to our first result, Theorem 1.2 
allows to consider also the case cr — t- oo which may be important for asymptotic considerations [21, 
20, 15]. 

Let us sketch the outline of the paper: We start with reformulating (l)-(2) as a fixed-point problem 
and then establish the unique solvability and a-priori estimates in L°° and L^ in Sections 3 and 
4, respectively. The proof of Theorem 1.1 is completed in Section 5 by extending these results to 
L^ 1 1 < p < cxD via interpolation arguments. In Section 6, we then construct a particular example 
showing that the estimate of Theorem 1.1 is sharp. Section 7 contains a short proof of Theorem 1.2. 
In Section 8, we derive estimates for derivatives and traces of the solution, and we introduce natural 
energy spaces for the problem (l)-(2). Finally, we present uniform estimates for the contraction factors 
of the source iteration which follow easily from the contraction estimates proved earlier. 

2. Reformulation as fixed-point equation 

Let us start by reformulating the radiative transfer problem as an equivalent integral equation in the 
usual way [10]. We define the scattering operator by 

/C(/)(r,v) :== / /c(r,v',v)0(r,v')dv', (r, v) e 7^ x V, (5) 

Jv 

further denote by 

(J5)(r_+tv,v)=e-J'X---+*^'")'^*g(r_,v), (r_,v)er_ (6) 

the extension of boundary values, and define a lifting 

£/(r_+tv,v)= / e-^>("--+'^^'")'^7(r-+sv,v)ds, (r_,v)Gr_, (7) 

Jo 

where < i < £(r-, v). By elementary calculations one can verify that 

(v • V + n)Jg = 0, (8) 

and 

(v • V + a)Cf = /, £/lr_ = 0. (9) 

This means that the extension Jg of the boundary values lies in the kernel of the differential operator 
and that the lifting £ is a right inverse of v • V + cr. The radiative transfer problem can then be seen 
to be equivalent to the following operator equation in integral form [10] 

(t> - CJC(t) + Cf + Jg. (10) 

To show the existence of a unique fixed-point, we will in the following sections select appropriate 
solution spaces, provide conditions on the data such that Cf and Jg lie in this space, and show that 
CK. is a contraction. 
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3. Solvability in L°° 

We will assume throughout that (A1)-(A3) hold and use the fact that for every point r £ TZ and any 
velocity v G V we can find a point (r_, v) on the inflow boundary r_ such that 

r = r_+tv with 0<i<^(r,v). (11) 

Also note that £{r,v) = i?(r_, v). We show flrst that CK, is a contraction on L°°{TZ x V). 
Lemma 3.1. For any (j) e L°°{Tl x V) there holds 

\\Cm\L^(n^v) < (l-e-""='"^~)||(/'|Uo.(Kxv). 
Proof. Using / = ICcf) in (7) and the assumption that a^ < a, we obtain for < t < £(r_, v) 

|(/:/C0)(r_+tv,v)|< / e- f^'+^^-^^^^^a'^ir^ + .sv,^r)ds Ml^^t^^^^ 
Jo 

Applying Banach's fixed-point theorem, we see that (10) has a unique solution (j> G L°°{TZ x V) 
whenever Cf and J'g are in L°°{TZ x V). This can be guaranteed by the following two results. 

Lemma 3.2. Assume that if G L°°(7^ x V). Then 

ll'C/IU°=(7^xV) < IK/llL-°(7^xV)• 

Proof. Using the definition of C, we obtain 

|A/(r-+iv,v)| < fe-f>^'^-+^^^-^''^r'\ef(r^ + sv,v)\ds<\\ef\\L^(nxV)- 

^' a 

Lemma 3.3. For any g G L°°(r_) there holds 

\\^I^Jg\\L°°(TZxV) < II^.9||l~(kxv) < ll5llL-(r_)- 

Proof. Since cr > we immediately obtain \jjg{r- + iv,v)| < |g(r_,v)|, which yields the second 
estimate. The first one follows from Lemma 3.1. D 

The proof reveals that g G L°°(r_) is in fact necessary to ensure that J'g is bounded. Combining the 
three previous Lemmas and the equivalence of the fixed-point equation (10) with the radiative transfer 
problem, we obtain 

Theorem 3.4. For any g G L°°(r_) and If G L°°{TZ x V), problem (l)-(2) has a unique solution 
(j) G L°°{TZ x V) which satisfies the a-priori bounds 

Uh^inxv) < ell^=^ll-=° (P/||l»(kxV) + ll5llLoo(r_)). 

Proof. The existence of a unique fixed-point for (10) follows from Lemma 3.1 and Banach's fixed point 
theorem. By the previous estimates, we get 

Ml^CizxV) < (1 - e-ll<'ll-~)||^||Loo(7^xv) + IK/IIl»(7^xV) + ll5llL=°(r_), 
from which the assertion is derived straight forward. D 

This completes the proof of Theorem 1.1 for the case p = oo. Note that actually no condition on the 
cross-section as was required here. 
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4. Solvability in L^ 

Setting w — v ■ V(f) + a(f) allows us to express the solution as (f) = Jg + Cw. The fixed-point problem 
(10) can then be stated equivalently as 

w = KCw + / + KJg, (j)^ Cw + Jg. (12) 

We want to show existence of a unique fixed-point for (12) in L^{TZ x V). To do so, we will first 
establish the contraction property for the operator iCC We will make use of the following well-known 
integral formula 

/ /(r,v)d(r,v) = / / /(r_+tv,v)|v-n|dtd(r_,v) (13) 

which directly follows from Fubini's theorem; see e.g. [32, 2, 14]. We can then show 
Lemma 4.1. For any w ^ L^{TZ x V) there holds 

WlCCw^.^nxV) < (1 - e-ll"=^ll-~)lklUi(7jxv)- 
Proof. By the definitions of /C and (jg , we get 

\\ICCw\\mnxV)< I I f k{r,v',^)\{Cw){r,V)\dv'd^dr 
JnJvJv 

cr,(r,v')|(/:w)(r,v')|dv'dr ==(*). 



inJv 
Using the definition of C and applying the integral formula (13) further yields 

(*)</" ( ' a,(r_+iv,v) /" e-J'>(''-+^^'^)^nw(i--+sv,v)|dsdi|n-v|d(r_,v) 
Jr^Jo Jo 

<[ f (l-e~ ■'■='*''"''" ''=('"-+''^^^)'*'')|w7(r_+sv,v)|ds|n-v|d(r_,v). 

Jr_ Jo 

Here we used cr^ < c and applied Fubini's theorem again to exchange the order of integrals with 
respect to ds and di and explicitly computed the latter. The assertion now follows directly. D 

A slightly weaker result was proven in a similar manner in [4]. To establish the existence of a fixed- 
point, we additionally have to require that / and ICJg are in L^{TZ x V). For the latter term, we 
use 

Lemma 4.2. For any g e i^(r_; |v ■ n|) there holds 

WICJgWL^nxV) < hsJgWmnxv) < (l-e-ll"='ll^~)||g||Li(r_;|v.„|)- 

Proof. By the definition of /C and ct^, we obtain 

\\ICJg\\mnxV)< I j (Ts{ry)\Jg{Yy)\dw'dv=\\a,Jg\\LnnxV)- 
JnJv 

Employing the definition of J" and the integral formula (13), yields 

hsJgWmnxV) < (l-e-ll"='ll^=°)||.9|Ui(r_;|v.„|), 

where in the last step, we used as < cf and a direct computation of the integral similar as in the proof 
of Lemma 4.2. D 

By Banach's fixed-point theorem and the previous estimates, we now obtain 
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Lemma 4.3. For any f G L^{TZ x V) and g G L^{r_; |v • n|), the fixed-point problem (12) has a 
unique solution w G L^{TZ x V) and there holds 

Mlhu.v) < e"^='"-=° (ll/liLl(7^xV) + (1 - e""'^='"^°°)ll.9llLi(r_;|v.„|))- 

To establish an L^ estimate for the solution (j) — Cw + Jg of problem (l)-(2), we have to establish 
additional bounds for Cw and J^g. 

Lemma 4.4. For any w G L^{TZ x V) and any g G L^{T^; |v • n|) there holds 

P~^/:w7||ii(7^xV) < \\w\\l^(1ZxV) and Wi^^Jgh^iUxV) < ll5llLi(r_;|v-n|)- 

Proof. These estimates follow from the integral formula (13) and direct computations. D 

A combination of the previous estimates now yields 

Theorem 4.5. For any f G L^{TZ x V) and g G i^(r_; |v • n|), the boundary value problem (l)-(2) 
has a unique solution (j) G L^{TZ x V) which satisfies 

IK^VIIlmt^xV) < ell"='ll^~ (ll/llLl(7^xV) + ll3llLi(r_;|v.n|)). 

Proof. The result follows from the representation (j) — Cw + J^g of the solution by applying the triangle 
inequality and using the estimates of Lemmas 4.3 and 4.4. D 

This completes the proof Theorem 1.1 for the case p = 1. Note that for our arguments, we did not 
use the condition on the scattering cross-section u'^ here. 

5. Solvability in L'p and proof of Theorem 1.1 

For establishing solvability in L^, we will utilize the results for L^ and L°° and the complex method 
of interpolation [7, 22] . Let us recall that for a cr- finite measure space (7W , d/i) 

P 

i.e., for 1 < p < oo the space L^ is an interpolation space between L^ and L°^ . In addition, the 
interpolation norm coincides with the norm of L^; see [22, Example 2.1.11]. 

We are now in the position to complete the 

Proof of Theorem 1.1. As a first step, let us establish the a-priori estimate for data that simulta- 
neously satisfy the requirements of Theorems 3.4 and 4.5. Noting that 

the a-priori bounds of these previous results can be written as 

\\<P\\LpCRy.V;e.-^) < e'^" (P/||Lp(7^xV;£-l) + ll.9llLp(r_ ;|v-n|)) for P e {1, oo}. 

Here Cp denotes the stability constant from Theorem 1.1. Using the linearity of the problem, we 
can decompose (f) = (j)g + (f)f, where (f)g and 0/ are the solutions of (l)-(2) with / = and g = 0, 
respectively. An application of the Riesz-Thorin theorem [7, 22] then yields 

I|0/I|lp(7?,xV;£-1) < e*^" IK.f I|lp(TCx V;i!-i) and \\(f)g\\LP(TZxV:e-^) < e"^" l|g||LP(r_:|v-n|) 

for any 1 < p < oo. From this estimate the a-priori estimate is derived via the triangle inequality The 
unique solvability for all admissible data follows by a density argument. D 
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6. Sharpness of the a-priori estimates 

In the following we show by example that the a-priori bound of Theorem 1.1 is sharp. Let TZ = .Bi(O) 
be the unit ball in M'' and let V = S^ he the unit sphere; note that hence v = v in the following. We 
consider the scattering operator 

}C(f){r, v) = cr(r, v)(/)(r, — v). 

This definition yields the essential property 

|/C(/)(r,v)| <o-(r,v)||</)||ioo(7^xv) 

An inspection of the previous results shows, that Lemma 3.1 and therefore all result of Section 3 hold 
true with a'^ — a also for this example. Equation (1) can then be written as 

V • V0(r, v) + cr(r, v) (0(r, v) - 0(r, -v)) = /(r, v). 

This construction yields that only directions v and — v are coupled in the transport equation. Let us 
fix one direction v and write 4>^{t) = 0(r- + ti'(r_, v)v, ±v). We then obtain 

- — </.+ + cr+(0+ -(/)") ==/+ in (0,1) with0+(O)=O, and 

-^^r +'T-(r-0+)-.r in (0,1) withr(l) = 0, 
and the solutions of these equations are given by 

(l)+{t) = </>-(l - i) = 1 - e-J'o'^^'^'^'^'. 
For I > 3 and k = 2'+'^, let us choose 

cr+(t) = cr"(l-f) = /^'(l-t)^ and 

f+it) = .f-{l-t) = a+(t)e-i'o"*'^+('')d^ 
After some basic calculations, one can see that 

||0^I|l~(o.i) = -/-^Cl) = no) = 1 - cxp ( - ■^) -: a{l), 
and that 

Wif^L^iUxV) = /+(0) ^ ./-(I) < cxp ( - (1 - l)2'('+3)-('+4)) 

= exp ( - (1 - y)|Ia±^||L"^^) =: ^(Oe-""*^". 
Note that a{l) and b{l) tend to one as I goes to infinity. Combining these estimates yields 

||0±||L=^=«(O>^ell-"^ll-||^/±|Uo.. 

This construction can be repeated for all directions v. Since we had a'^ = a here, this shows that the 
estimate of Theorem 1.1 is sharp at least in the case p = oo. 

7. Proof of Theorem 1.2 

We now illustrate that better a-priori estimates can be obtained, if some absorption is present. Let 
the assumptions of Theorem 1.2 hold. Formally multiplying (1) with (/)|(/)|p^^ yields 

P 
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By integrating this equation over TZ x V, performing integration by parts, using the boundary condi- 
tions (2), and rearranging terms, we obtain 

1 

P 



lk'''^llLP(TCxV) + -|l<^llLP(r+;|v.„|) 



, ,^-2/C0d(r,v)+ / fcf,\<l>r'd{r,v) + -\\gr 
The first term on the right hand side can be estimated with Holder's inequality by 

l-RxV 

For the last step, we used the following basic estimates for the scattering operator 

\\IC<I>\\l^<\\-\\l^\W4>\\l^ and \\ilC<j>\\L^ < W^U^Uh^ 
a a a 

from which one obtains by interpolation that \\a ? /C0||lp < ||cp0||lp. To bound the term involving 
the right hand side /, we apply Holder's and Young's inequality, to get 



TZxV 



<r.^-pi\\a'^f\\l^^^^^^+.P-^\\ah\\Unxvr 



Putting all estimates together and multiplying by p, finally leads to 

''IW^^WIhUxV) + ll'^liLP(r+;|v.n|) ^ ^'^'^ 1^^ /H Lp(TCx V) + ll.9liLP(r_ :|v.„|) ' 

From this estimate, the assertion of Theorem 1.2 now follows directly. 

8. Additional results 

To complete our discussion, we collect in the following some further results which follow more or less 
directly from our previous considerations. 

8.1. Estimates for the derivatives. Using the a-priori estimates of Theorem 1.1 and the fixed-point 
equation (10), it is straight-forward to obtain also estimates for the directional derivatives v • Vi/). Let 
us first consider the case p ~ 1, where we have 

Lemma 8.1. Under the assumptions of Theorem 4-5 one has 

l|v • V^IUnKxV) < 2ell-=^ll-- iWfWmnxV) + ll.9llLi(r_;|v.„|))- 
Proof. By (j) = Cw + Jg and the properties of the operators C and J , we obtain 

l|v • V(/)||li(7j,xv) < \\w\\L^(nxV) + IIo-^IIlicrxv)- 
The first term can be estimated by Lemma 4.3, and for the second, we use 

\W<P\\L^(nxV) < \\<^^w\\L^nxV) + IWJghHTZxV) < II^IIliCRxV) + l!.9l|Li(r_;[v-n|)- 

The second estimate for the boundary term is obtained as in Lemma 4.2. D 

For the case p — oo, we have 

Lemma 8.2. Under the assumptions of Theorem 3.4 there holds 

pv • V^L^inxv) < (1 + 2|la^|U^)ell^=^ll-~(-xv) {\\fh^(^nxv) + lig||L-(r_;iv.„i))- 
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Proof. The identity v • \7(f> ~ ICcp — acf) + / yields 

pV • V0||loo(kxV) < (IKc^!,||L°o(7^xV) + IKcr||Loo(Kx V) ) II0I|L~(7^X V) + IK/I|l~(Kx V) ■ 

The estimate then follows from the bounds of Theorem 3.4 and the condition a'^ < a. D 

Arguing as in the proof of Theorem 1.1, the case 1 < p < oo is then covered by 
Theorem 8.3. Under the assumptions of Theorem 1.1 there holds 

pi-H • V4>\\LHn^v) < 2(1 + lk£||L»)e^-(pi-t/|Up(^xV) + llffllLp(r_;|v.„|)). 
Again, we can obtain in a similar way stronger estimates under additional assumptions. 
Theorem 8.4. Let the conditions of Theorem 1.2 hold. Then 

\\a-v-^w ■ V</>|Up(kxV) < (2j^"' + l)lk^"Vl|Lp(KxV) + 2i^^/P||5l|Lp(r_;|v.„|))- 
The use of weighted norms again substantially simplifies the derivation of these results. 

8.2. Energy space and a trace lemma. The norms in which we obtained the a-priori estimates of 
Theorem 1.1 and Theorem 8.3 suggest to define the following energy space 

WP = {0 : 7^ X V ^ R : ^"^0 G 1^(11 x V), t^^<r ■ Vcj) e LP{TZ x V)}. (14) 

The natural norm for this space is given by 

For functions in the space W^, wc have the following result for traces. 

Theorem 8.5. The trace operators "f± : W^ — > LP{T±: |v • n|) are continuous and surjective. 

A proof of this statement follows easily by direct computation; see also [8]. For other estimates and 
general material about traces for radiative transfer problems let us refer to [13, 23, 14]. 

Using the results of Sections 3-5, we also obtain that 

\\e'-^ {^-V^ + a^,- /C0) ||LP(7^xv) < (1 + 2||a^|U^)||</-||w- 
The individual operators could be estimated in the same way. Summarizing, we obtain 
Theorem 8.6. Let (A1)-(A3) hold. Then the mapping 

WP ^LP{nxV-JP-^)xLP{T_-\w-n\), (P ^ {v ■ Vcj) + cr4> - ICc/), -/^c/)) 
is continuous and boundedly invertible. 

This result shows that the assumptions on the data cannot be relaxed when searching for solutions 
in the energy space W^. Under the stronger assumptions of Theorem 1.2, we can define in a similar 
manner an energy space 

WP = {crp(/) e LP{TZ X V), cr^-^v • V(f> e LP{n x V)}. 

Results analogous to Theorem 8.5 and 8.6 can easily be derived also for this space. For the correspond- 
ing statements it suffices to replace the weight function £ by cr~^; compare also with Theorem 1.1 and 
1.2. 
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8.3. Spectral estimates and convergence of the fixed-point iterations. The solvability results 
of the previous sections were based on Banach's fixed-point theorem. The corresponding fixed-point 
iteration reads 



l^n+l 



CKcj),, + Cf + Jg. (15) 



We show now that under our general assumptions (Al)-(A3), the spectral radius of the fixed-point 
operator CIC is always uniformly bounded away from one. 



Theorem 8.7. Let (A1)-(A3) hold. Then for all 1 < p < oo 



PpiCJCj := lim x/WiCICy'WLPiKxv-.e- 

n— 7-00 V 



< 1 - e^^" 



Proof. The case p = oo follows immediately from Lemma 3.1. For p — 1, on the other hand, we can 
estimate the powers of the fixed-point operator by 

||(£/C)"||l1(7JxV;£-1) = P"^('^^)"^IIl1(7?,x V) < IK^^-^IIlMT^-X V) I1^-CHli^(7J,x V) li^^lli'(Kx V) • 

The first two terms can be bounded by Lemma 4.2 and 4.1, and for the third term we use the estimate 
II^^IIli(TCxV) ^ II'''s^IIl°°(TCxV)- From this we obtain the estimate for the spectral radius for p — 1. 
The general case then follows again by interpolation arguments. D 

Our analysis thus shows that under the weak sub-criticality assumptions (A3), the source iteration 
(15) converges in L^ for any 1 < p < oo with a contraction factor 1 — e^'-^p. Note that no positive 
lower bounds on the absorption are needed for the convergence. The same arguments may be used to 
analyze other fixed-point iterations, cf. [1]. 
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